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ABSTRACT 



We study the final stages of the common envelope (CE) evolution and find that a substantial fraction 
of the ejected mass does not reach the escape velocity. To reach this conclusion we use a self-similar 
solution under simplifying assumptions. Most of the gravitational energy of a companion white dwarf 
(WD) is released in the envelope of a massive asymptotic giant branch (AGB) or the red giant branch 
(RGB) star in a very short time. This rapid energy release forms a blast wave in the envelope. We 
follow the blast wave propagation from the center of the AGB outwards, and show that ~ 1-10 per 
cent of the ejected envelope remains bound to the remnant binary system. We suggest that due to 
angular momentum conservation and further interaction with the binary system, the fallback material 
forms a circumbinary disk around the post- AGB Core and the companion WD. The interaction of the 
circumbinary disk with the binary system will reduce the orbital separation much more than expected 
from the dynamical phase (where the envelope is ejected) of the CE alone. The smaller orbital separation 
favors a merger at the end of the CE phase or a short time after, while the core is still hot. This is 
another channel to form a massive WD with super- Chandrasekhar mass that might explode as a type la 
supernova. We term this the core-degenerate (CD) scenario. 

Subject headings: stars: AGB and post-AGB - (stars:) binaries: close - (stars:) supernovae: general - 
(stars:) white dwarfs 



INTRODUCTION 



The 35 years old ( Paczvnski 19761) common envelope (CE) model is in the heart of the formation of many close 



binary systems (e.g. 



2001; Wcbbink 2008 



Mever fe Mever-Hofmeiste3 ll979i llben fc Liviolll993t iTaam fe Sandauistllioool IPodsiadlowski 



Taam fe Rickerll2010l ). The CE is defined as a structure where two stars share the envelope. 



In most cases the mass of the envelope comes from a giant star as it engulfs its companion, which is either a main- 
sequence (MS) star or a white d warf (WD). During t he CE phase the orbital separation decreases due to gravitational 
drag and tidal interaction (e.g.. llben fe Liviolll993f) . The transfer of orbital energy and angular momentum to the 
envelope, as well as other possible energy sources, lead to the ejection of the envelope. 

One of the major unsolved questions of the CE phase is the final orbital separation. It is customary to equate 
the gravitational energy released by the spiraling-in binary system Eq, to the binding energy of the envelope -E^ind- 
An efficiency parameter «ce is introduced such that the final orbital separation is derived fr om the equality -Ehinri = 
oiceE q. However, different studies deduced different values of acE for different systems. Ivanova fc Chaichenetsl 
(|201l[ ) suggest that the enthalpy rather than the internal energy should be included in calculating the binding energy. 
A better understanding of the physical processes that determined the final orbital separation is required. 

Most of the energy is deposited into the env elope in a very short time, as the spiralin g-in process is accelerated 
with the decrease of the orbital separation (e.g., Rasio fc Livio 1996 : Livio fc Soker 19881) . There are also allegedly 
contradicting re s ults, w hich suggest that the time for energy release in the envelope is relatively long. For example, 
Sandquist et al.1 (|l998l ) found that the energy is released mainly during a stage lasting ~ 200 days (see their fig. 
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11). iDe Marco et al.l (|2003l ) found that in representative cases the entire timescale for CE evolution is 9-18 years. 
However this timescal e was obtained by waiting for a negligible amount of material to remain in the envelope 
(|De Marco et al.ll2009l ). and therefore it represents a much longer timescale than the strong energy release timescale. 
This may also suggest that the s trong energy release timescale is of the order of ~ 1 year or less. Recent numerical 
simulations of iPassv et al.1 (|201lf ) also show that the strong energy release occurs in a timescale of 100-200 days. 



The energy that is released during the CE is usually composed of more than one bur st, and it takes for th e 
companion a few orbits to strip the asymptotic giant branch (AGB) star envelope away (e.g., Sandauist et al. 19981) . 
If most of the energy is released within a short time deep in the envelope, i.e., shorter than the dynamical time 
of the outer parts of the envelope, then one can speak of a burst of energy. Such a burst is similar in some of its 
properties to a blast wave which propagates outwards to the stellar surface. This is the scenario we discuss in the 
paper. 



Sandauist et al. ( 19981) showed that some of the mass of the envelope can remain bound to one or both of 
the interacting stars. They simulated a 5 M Q AGB interacting with a 0.6 M ff companio n, and found that the 
companion unbinds ~ 1.55 M Q (~ 23 per cent) of the AGB envelope. IPassv et al.l ([201 ll ) present more extreme 
results in their simulations. They find that when the envelope is lifted away from the binary, > 80 per cent of the 
envelope remains bound to the binary. However, in their simulations rotation of the AGB was not included, so 
the material that remains bound is probably overestimated. They conclude th at in some cases parts of the AGB 
envelope remains bound or marginally bound to the remnant post-AGB core. ISandquist et al.l ([19981 ) also found 
that a differentially rotating structure resembling a thick disk surrounds the remnant binary during an intermediate 
phase of the CE interacti on. This stag e occurs briefly befor e energy deposited by the companion and remnant AGB 
core drive the m a ss aw ay. ISokerl rtl992j ) (see also lSokerEool have analytically obtained a similar thick disk structure. 
De Marco et al.l ([201 11 ) suggest that envelope material which is still bound to the binary system at the end of the 
CE, will fall back onto the system and will form a circumbinary disk. They suggest that such a disk might have 
som e dynamical effects on t he binary period. Other simulations show a post-CE gas to be bound to the companion, 



e-g- ; 



Lombardi etaH |2006) 



In this paper we discuss the final stages of the CE under the assumption of a rapid binary-gravitational energy 
release, and the formation of a circumbinary disk. We suggest that such a scenario can result in an early merger of 
the WD with the post-AGB core, that will explode as a Type la supernovae (SNe la) - the core-degenerate (CD) 
scenario. 



Sparks fc Stecherl (11974 ) develope d a me rger scenario of a WD with an AGB core, in which the result is a type 
II SN. In a later study, iLivio fc Riessl (|2003l ) found that the merger of a WD with an AGB core leads to a SN la 
that occurs at the end of the CE phase or shortly after. An important conclusion of Livio fc Riess ( 20031) is that for 
merger to occur the AGB star should be massive. In the present paper we also reach the same conclusion. 

We discuss the problem of the energy budget and characteristic timescale at the end of the CE in section [2] In 
section [3] we discuss the properties of the blast wave. For such a disk to form some of the material needs to fall-back 
onto the star. We follow the blast wave and show that it indeed leaves some of the envelope material bound to 
the binary system. We discuss the formation of a circumbinary disk and its potential consequences on the binary 
system and its final fate. We argue that in many cases such a disk induces a merger and discuss the implications of 
the formation of such an early merger (section |4|) . The most interesting implication is the CD scenario we propose 
for a SNe la. In the CD scenario a merger of a WD companion with the post-AGB core occurs while the core is still 
hot, and later explodes as a SN la. In other cases the merger might collapse to a neutron star. We summarize in 
section [5] 
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2. PROBLEM SETUP 

2.1. Energy considerations 

In this section we discuss the energy considerations behind the process of expelling the envelope. We follow 
a stellar companion spiraling-in inside the AGB envelope, and concentrate on the orbital separation a/ where the 
liberated orbital gravitational energy is about equal to the envelope binding energy. As stated above, the final 
stages of the spiraling-in occur fast. Even if the envelope inward to the companion reaches synchronization with 
the orbital motion, the spiraling-in will continue. The reason is that the envelope outer to the companion cannot 
be synchronized, and tidal interaction are very strong. This tidal interaction is similar, but not identical, to that 
with a circumbinary disk. As we will see in section |4l such interactions are very efficient in reducing the orbital 
separation. Based on this discussion, we assume that an energy about equal to the binding energy of the envelope 
is released in a very short time. We examine two analytical prescriptions for this energy deposition in section [3] 

We start by examining the involved energies. For our typ ic al giant structur e we a pproxi mate t h e AG B envelope 



structures from the models of iNordhaus fe Blackmarj ([20061 ) . iTauris fc Dewi (l200lh . and ISokerl ([19921 ). that are 



qualitatively similar. As we are aiming a t performing self similar calculations, we approximate the envelope density 
profile by a power law, e.g.. lSokerl (|1992l ) 

with uj = 2. The mass inward to radius r is given by 

M (r) ~ M corc + M cnv (r) = M corc + f ^r 2 p dr. (2) 

Our AGB stellar mass is M* = 5 M,®, its core mass is M corc = 0.77 M,®, its envelope mass is M env = 4.23 Mq 
and the stellar radius is i?+ = 310 Kq. For these parameters and oj — 2 we find from equations ([T|) and @ that 
A ~ 3.1 x 10 19 g cm -1 . The core radius is very small, fewxO.Ol R©. 

The binding energy of the envelope mass outside radius r is given by 

[ R * G(M(r) + M 2 ) 2 , 
E hind = / 1 1 '- 2 Mnr 2 pdr, (3) 

J T ' 

where we have taken into account the companion mass M 2 , and neglected thermal energy of the gas. Before the 
companion enters the envelope its influence is mainly through tidal interaction that makes the envelope spin-up. 
This spin-up and the entrance of the companion to the envelope is likely to cause some mass to be lost. This 
mass will be signific ant if the companion manages to bring the envelope to synchronization with the orbital motion 
dBear fc Sokerll2010h . We consider a case where if a synchronization is achieved it does so before the AGB reaches 



its maximum radius. As the AGB continues to expand it swallows the companion. We assume that not much mass 
was lost before the onset of the CE. The gravitational energy released by the companion inside the envelope as it 
spirals in from the stellar surface to binary separation a is given by 

GM(a)M 2 GM*M 2 

EG = 2a 2RT- (4) 

The binding energy and the released gravitational energy for our AGB model and for M 2 = 0.6 M Q are plotted in 
Fig. [TJ For an ~ 1.4 R the two energies become equal to each other and to Ef ~ 6 x 10 47 erg. If we use the 
alpha prescription, acE-E-Gfafa) = -Ebind^fa), then the ejection of the envelope occurs at an envelope separation of 
af a — acEafi- 



2.2. Timescales 



In this section we examine some timescales relevant to the final stages of the common-envelope phase. The 
spiraling-in timescale due to tidal interaction when the companion is inside the envelope is about equal to the 
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Fig. 1. — The binding energy of the AGB envelope mass residing at r > a, and the released gravitational energy, for our case of a 
5 Mq AGB with envelope density ~ r~ 2 and a 0.6 Mq companion, as function of the orbital separation. As the companion spirals-in 
towards the center it releases more and more energy, namely is releases most of its gravitational energy close to the center. Until it 
reaches the inner intersection of the binding energy of the AGB and the released gravitational energy (af ~ 1.4 Rq), the companion 
liberates Ef ~ 6 X 10 47 erg and unbinds most of the AGB envelope. 



circularization time 



rtid(a) 



z (a) ~ 6.9 



2 x 10 4 L 



300 R Q 

M 2 
0.2M(a) 



/ -Mcnv(a) 

V 4M 
-l 

' 1 



M 2 



M(a) 



M(a) 

5 M 

1 



(5) 



where L* is the stellar luminosity, and the circularization time is from IVerbunt fc Phinnevl (| 1995! ) . and where inside 
the envelope we take the radius inward to the companion equals the orbital separation r = a. In this expression only 
the tidal interaction with the envelope inward to the orbit is considered, however there is also a tidal interaction 
with the envelope exterior to the orbit (see section POT) . that shortens the process. 

We assume that the convection in the envelope inner to the companion radius still exists. For that, the tidal 
dissipation is assumed to be valid. Spiral wave dissipation exists as well (e.g., Sandauist et al. 19981) . for which the 
characteristic timescale is ~ 200 days. 

There are two additional timescales to consider. These are the Keplerian timescale 

„3 



T K {a) = 2n 

and the friction timescale (e.g., Iben fc Livio 19931) 

Tfric(a) = 

where 



GM(a 
GM{a)M 2 

2aLdrag(o) 



idrag(a) = ^Lc^K 

is the luminosity due to energy loss by drag, v^(a) = {GM(a) / a) 1 / 2 is the Keplerian velocity, 

2GM2 

-ftaccla] — t 



(6) 

(7) 
(8) 

(9) 
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Fig. 2. — The characteristic timescales for unbinding the AGB envelope, for our case of a 5 Mq AGB with envelope density ~ r~ 2 
and a 0.6 companion. T t ^: the spiralling-in timescale due to tidal interaction with the envelope residing inward to the companion 
orbit. Tf r ; c : the friction timescale. tjj : The Keplerian timescale. r: the relevant timescale at radius r according to equation fl 1 1) . 



is the Bondi-Hoyle accretion radius, and we take the constant £ = 1. The usage of i>k in the equation (J5|) is 
an approximation. A more accurate expression takes into account the rotation of the envelope and the sound 
speed. Including envelope rotation reduces the relative velocity, while the sound speed increases it. Overall, the 
approximation used here is adequate. 

The accretion rate might be less than the Bondi-Hoyle-Lyttleton value. However, in the drag equation there is 
another term resulting from interaction of the companion with mass that is not accreted. This term has the form 
of (H|), but with the coefficient ln(_R max /i ? a cc) instead of £, where i?m ax is the distance up to where the companion 
gravity influences the envelope gas (e.g., Ruderman fc Spiegell 1971 ). This implies that even if the accretion rate 



is much below the Bondi-Hoyle-Lyttleton value, the term we ignored will become important and the value of the 
drag will be about the same as in equation ([8]) with £ ~ 1. Therefore, even if there is a large deviation from 
Bondi-Hoyle-Lyttleton accretion rate, the value of rf r j c will only change by a small factor from the value used here. 

The spiraling-in timescale at orbital separation a is dictated by the most efficient process, but under our 
assumptions it cannot be shorter than the Keplerian timescale. It is given by 

r(a) = max[min(r t id, r fric ), r K \. (10) 

The value of r(a) is plotted in Fig. [2 for our stellar model and for M2 = 0.6 M . We note that the tidal interaction 
is small compared with the gravitation drag inside the envelope. At the radius where the released gravitational 
energy equals the envelope binding energy, an, t? = r(afi) ~ 0.95 days. From Fig. [I] it is evident that most of the 
gravitational energy is released very close to an. For example, 90 per cent (80 per cent) of the energy is released 
when the companion spirals-in from 14 R (7.1 R©) to Of, and the timescale in that location is t(14 Rq) ~ 6.1 days 
(t(7.1 R q ) ~ 2.4 days). Overall, most of the energy is deposited in the inner region of the envelope and during a 
time much shorter than the dynamical time of the outer regions. The consequences of this instantaneous energy 
liberation are discussed in the following sections. 

We note that had we taken the spiral wave dissipation timescale instead of the tidal timescale (e.g., from fig. 3 



Sandauist et al.l (|1998T )). we would have gotten a shorter timescale than Tya (the blue line in Fig. [2] would have 



been located a little lower). But, according to equation (fT0|) the relevant timescale for r would not have changed. 
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Fig. 3. — The asymptotic self-similar solution of an intense explosion blast wave which satisfies u = (7 — t)/(7 + 1): as solved bv lSedovl 
j 19591) . 

3. THE BLAST WAVE PROPAGATION THROUGH THE AGB ENVELOPE 

We perform a self similar calc ulation of a bl ast wave propagating from the interior of the AGB outwards. The 
calculation is based on the work of lSedovl dl959h . To facilitate the analytic solution we consider two limiting cases 
of the short energy deposition discussed above. In the first case we assume an instantaneous release of the binary 
gravitational energy (section 13. 1|) , while in the second case the energy is released at a constant rate over a short 
period (section l3~!2|) . 



3.1. Instantaneous energy release 

H ere w e assume that the energy is deposited instantaneously at t = 0. In this case (based on chapter 14.4 in 



Sedovlll959l ) we neglect the initial pressure in the AGB envelope and consider an adiabatic motion of a perfect gas. 
The density profile before the blast is given by equation(fT]), and the envelope is static. The general solution for the 
blast wave is rather complicated for a general value of u), and adiabatic exponent 7. However when the relation 
uj = (7 — + 1) holds, as in the case of the AGB model we use with oj = 2 and 7 = 5/3, there is a simple 

asymptotic solution for the velocity, density, and pressure behind the shock, given respectively by 

x 2 r 

V = D2A 



37- It 



Ah + 1) a 



where the index 2 refers to the immediate postshock values, and where A is the self-similar variable 



A =U) rt "-s- (12) 

with T2 the location of the Shockwave, and Eq the energy deposited in the envelope. The constant a is determined 
from energy conservation. The self-similar solution for u> = 2 and 7 = 5/3 is shown in Fig. [3] 

In addition to the energy deposited by the spiraling-in companion, the envelope itself has thermal energy i?th, 
which we also add to that energy. We use the virial theorem E t h = 0.5i?bmd to find the energy deposited in the 
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cnvclope 

3 

Eo = -Ebind + E t h — --Bbind (13) 



The justification for considering the thermal energy of the envelope was recently discussed bv lDe Marco et al.l (|201lh . 
We use this energy in equation (fT2]) to obtain the solution that is shown in Fig. 01 

We use the solution to estimate the mass that stays bound to the remnant binary system. We calculate the 
escape velocity Uesc(^) = [2G(M(r) + M2)/?-] 1 / 2 taking the companion mass into account, as its mass contributes to 
the gravitational force at the time the enve lope is ejected. Radiation pressure is not an important factor in expelling 



material outwards (e.g., iPassv et al.ll201l[ ). We therefore neglect it, and set the condition for ejecting material to 
have its velocity larger than the escape velocity. The amount of material that does not reach the escape velocity and 
hence remains bound is determined from the condition at the end of the calculation v(r, tend) < v esc (r,t en d), where 
t cn d is the time when the Shockwave reaches i?*. This material is expected to fall-back towards the center, as the 
pressure decreases inwards (Fig. |4|). For our case we find i cn d = 9.5 days ~ 10rf. This shows that our assumption 
that the energy is released instantaneously is quite adequate here. 

We find that for = 0.6 M Q , our characteristic case, the fall-back mass is Mfb ~ 0.22 M Q . For larger M2 
more of the envelope mass is falling back, as we show in Fig. [5j In Fig. [5]we also show the ratio between the time it 
takes for the blast wave to reach the AGB radius (i e nd) and the spiraling-in timescale at the final orbital separation 

(Tf). 

The formation of accretion disks around one or two of the post-CE stars by fall back m aterial in the post- AGB 



phase was discussed before, but at much later times after the end of the CE phase (e.g., ISokerl 120011) . Here we 
discuss a formation of a disk immediately as the main CE phase ends, and emphasize a disk around both stars-a 
circumbinary disk. 

The total internal energy of the gas should include ionization energy as well, as the gas cools and recombines 
when it is expelled. Here we include only the thermal energy. The ionization energy might b eco me an important 



sourc e in extended AGB stars that have low binding energy per unit mass fe.g.. lHan et al.ll2002t see lSoker fc Harpaz 



2003 for a different view). For AGB stars the effects of ionization energy in the internal energy of the gas may 
be important. As we consider here massive AGB stars, it seems the ionization energy will not change much the 
amount of expelled mass. While the binding energy goes as -Ebind oc M 2 , the ionization energy goes as E\ oc M. 
For an envelope of M onv = 4.23 M Q the ionization energy available is ~ 10 47 erg, while the energy deposited in our 
calculation is ~ 9 x 10 47 erg. Furthermore, the bound mass resides in the inner part of the expanding envelope. The 
recombination radiation of the bound mass will help in part to expel the mass further out. Therefore, the ~ 11 per 
cent additional energy of recombination will not be efficiently deposited into the bound mass itself. It might expel 
an extra small amount of mass that is just barely bound. 

In our analytic study we are limited to spherically symmetric geometry. In asymmetrical cases we expect that 
more mass will stay bound. The reasoning goes as follows. Part of the mass stays bound because a fraction of the 
envelope leaves with much more energy than is required to escape, hence leaving part of the envelope with a total 
negative energy. An asymmetric effect, either a denser envelope in the equator or more energy injection into the 
equatorial plane, for example, will increase the uneven energy distribution in the envelope. More energy deposition 
in the equatorial plane might leave more mass bound near the polar directions. The question is what is exactly the 
value of the specific angular momentum of the bound envelope. It is possible that interaction of the bound mass 
with the binary will lead to the formation of a circumbinary disk even when the initial specific angular momentum 
is small. 
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3.2. Continuous energy deposition 



We find that the timescale for binary energy deposition rp is smaller but not negligible relatively to the time for 
the shock to reach the stellar surface tend- For that we now turn to calculate a case where the energy is deposited 
continuously during a finite time. In this second case we do not neglect the press ure of the eny elope The self similar 
calculation of the blast wave we perform in this section is based on chapter 6.6 of lSedovl (|1959f ). For this self similar 
solution the energy deposition rate is Eq oc t 2 ( 5 ~ 2 ")/ w oc t, where the last equality is for oj = 2. 

As a result of the continues energy deposition and the counter pressure of the pre-shock envelope, a void is 
created in the inner region, and the post-shock material is contained within the blast shell, defined as the volume 
between the outer boundary of the void and the Shockwave. The blast shell moves outwards, until it reaches the 
stellar surface. The problem with this description is that in our case the energy source is the companion. When the 
inner boundary of blast shell crosses the companion, there is no mechanism to transfer the energy from the central 
engine to the shell, because the companion is inside the void. We therefore limit ourselves for treating a blast with 
a continuous energy release to low radii and short times, when the blast shell is close to the center. As in our case 
the continuous energy release solution is relevant only to short times, we stop the energy repositioning at Tf . 



The self similar variable here is 



A = r(pAG)^t^ = —. 

T2 



(14) 



The self-similar functions are R, V, Z and A4, defined as (chapter 6.6 in lSedovlll959h 



P{r,t) = 




v(r,t) = 


r 


p(r,t) = 




M(r,t) = 


7& M ™ 



(15) 



jGt 4 ' 



The differential equations that describe the blast, in terms of the self similar functions are 
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and the boundary conditions at the shock front are 

2g( 7 + l) (w-1) (w-3) 



i?2 

^2 = - 



7T cj 2 7 (2 + 7 — 1) 
4(9-1) 



-(7 + 1) 

4(g( 7 -l)-2 7 ) (2g + 7 -l) 1 J 

Z>9 = 



(7+I) 



A4 2 = 

where 



8q (lu - 1) 



w 2 7 



7TW 2 7 

g= ~2(^-l) (w-3)/3' (18) 
and /3 is a constant determined from energy conservation. 

The solution is plotted in Fig. [5] We can see that the blast shell resides in the inner part of the envelope, and 
supersonically moves outward. The energy in the blast wave is the same as that in the previous case (equation 1 13|>. 
such that we are back to the intense explosion case (section [3. ip . Namely, if the process does not last for a long time 
(see below), the two cases studied in this section lead to the same result that a non- negligible mass stays bound. 



POSSIBLE IMPLICATIONS 



4.1. The role of a circumbinary disk 

The fall-back mass is expected to posses angular momentum, originally coming from the orbital angular mo- 
mentum of the companion. We suggest that the fall-back mass forms a disk around the remnant binary system 
- the companion and the leftover core of the AGB. Even if the specific angular momentum of the fall back mass 
is not sufficient to form a circumbinary disk, later interaction of this gas with the binary system might lead to 
the formation of a disk. Such an interaction is thought to lead to the formation of the c i rcumbinary disks around 



post- AGB binary stars w ith an orbital sepa r ation of ~ 1 AU ([van Winckel et all bood Ivan Winckel et al 



2006 



van Winckel et al.l 120091 ) . Ivan Winckel et all (l2009h conclude that circumbinary disc can tremendously affect the 
evolution of the binary system in the center, and that the formation of circumbinary gravitationally bound disks 
occur in a wide range of post-AGB binaries. In these binaries the typical orbital separation is ~ 1 AU, many of 
them have substantial eccentricity, and even if they avoided a CE phase, the companion had a strong interaction 
with the AGB envelope. 



In the simulations of iPassv et al.l (| 2 llh even larger larger than found in our calculations amount of material 
remains bound to the binary system. Howe ver, they find that the orbital separations of the post-spiral-in phase are 
larger than the observed ones. Though that IPassv et al.l (|201ll ) found that there is enough orbital energy to unbind 
the envelope, the binary system in their simulations remains in a completely evacuated volume and the leftover 
bound envelope is at ~ 100 R Q . 

We therefore assume that the fall back gas forms a circumbinary disk. The binary system transfers angular 
momentum to the circumbinary dis k, and consequently the r adius of the disk expands, the binary separation decreases 
and the eccentricity increases (e.g.. lArtvmowicz et al.lll99lh . We assu me that the structure of t he circumbinary disk 
for med in our scenar i o rese mbles the structure of the disk studied in lArtvmowicz et al] (|199lf ). In the calculations 
of lArtvmowicz et al] (|199ll ) the disk extends from the nearest stable circumbinary orbit of ~ 2.5a up to 6a, where 
a is the binary separation. We substitute for the binary system used in section [3l and take the binary system to 
be composed of M core = 0.77 M Q and M-i = 0.6 M components at a separation of a,{ ~ 1.4 R©, which is the final 
binary separation of the CE phase. We assume that the initial binary orbit is circular (eo = 0). The circumbinary 



- 10 - 



di sk mass is the mos t unce rtain quantity here. We here scale it with A/disk, o = Mfb — 0.2 M Q . Scaling the results 
of lArtvmowicz et al.l (|199lh by these values we find that the rate of change of the semi-major axis is 

(<£) (s^FTizr) < 19 > 



a 

- ~ -3 x 10" 
a 



and the rate of change of the eccentricity is 
where 

n b = ^ Wog+ MsT ^ 4 2 x 10 _ 4 s -i (21) 

is the binary orbital angular velocity, and 

qd - M corc + M 2 (22) 
is the disk to binary mass ratio. At the end of the CE, when the circumbinary disk -binary interaction starts, the 
disk to binary mass ratio is g^.o — 0.16. 

Another effect we take into account is the dissipation of the circumbinary disk. The energy released by the 
shrinking binary orbit is transferred to the disk. Because the total power is larger than the Eddington limit, a large 
fraction of the energy is carried by gas escaping from the disk. We take this mass loss into account by equating 
the energy carried by the wind to that released by the shrinking binary system. As material probably escape from 
the entire disk, we take as an approximation the wind to leave from the middle of the disk, i? w i n d(i) = 4.25a(i). 
More likely, more mass will leave from the inner region of the disk, hence less mass can carry the same amount of 
energy. This will leave more mass in the disk, making the scenario discussed here more efficient even. In general, 
winds escape from stars and disks with a velocity about equal to the escape velocity. Therefore, the total energy 
carried by the wind is the sum of the binding energy of the gas and the kinetic energy at infinity. We consider also 
the kinetic energy of the wind, assuming its terminal velocity is the escape velocity from i? W ind- The kinetic energy 
introduces a factors of \ in equation ((26)) below. 

We use equations (fT9|) and ([20)) to integrate the semi-major axis decrease and eccentricity increase due to the 
circumbinary disk - binary interaction. We take the circumbinary disk mass to reduce according to the released 
orbital energy from the reduction of the binary orbit 

GM coie M 2 GM corc M 2 
dE{t) = + 2a(t) ■ (23) 

The binding energy of the escaping disk material per unit mass is 

<iE'bind,disk(£) G(M core + M 2 ) 



dm 2i? w i n( j 

and the kinetic energy is parameterized as 



(24) 



dm dm 

Hence, the rate of change of the disk mass is 



^-^kin.disk / -, \ ^-^bind,disk /r»r-\ 

1 • (25) 



X dm 



Our results for an initial eccentricity eo = are presented in Fig. [7] The decrease in orbital separation and 
increase in eccentricity are clearly seen. The circumbinary disk - binary interaction terminates when one of two 
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things occur: either the circumbinary disk is completely depleted, or the two stars merge (occurs when e ~ 1). We 
term the time period of the circumbinary disk - binary interaction T\, As can be understood from Fig. [7J taking 
eo > would shorten t\ , so taking eo = is a conservative parameter selection. As stated above, the most uncertain 
parameter is the initial disk mass, which sets the value of qd,o (equation |22|). We repeat the calculation for a few 
values of qd,o- For the case where the disk wind is ejected at the escape velocity (namely x — 2) we find that the 
circumbinary disk - binary interaction leads to a merger for g^.o ^ 0.12, without the necessity of a following process 
(emission of gravitational waves (GW) , discussed in section 14.31 below) . 

We note that post-CE circumbinary disks were suggested u nder other as sumptions about the CE evolution. The 
formation of a post-CE thick circumbinary disk was invoked bv lSokerl |1992) in a case of a more gradual spiraling-in 
process where part of the envelope remains in an hydrostatic equilibrium rather than falling back. A differentially 
ro tating thick disk or tor us was found in intermediate stages of the CE evolution in the 3D numerical simulations 
of lSandauist et al.l (|1998[ ). In these cases a rapid merging is expected as well. 



Iva nova ( 20111) suggests another mechanism for post-CE merger, applicable for stars with non-degenerate cores. 



Ivanoval (|201lf ) divides the CE phase into two phases. In the first stage the envelope is expelled, and a remnant 
is left in the center. This stage takes ~ 1 yr. Some envelope mass is left above the core. In the second stage the 
left-over envelope under goes thermal readjustment, and its outer region expands. This occurs on a timescale of a few 
x 10 3 yr. Ilvanoval (|201ll ) concludes that in many cases the left-over envelope would expand beyond the companion, 
forming another CE phase that substantially increase the merger probability. 

As we found above that some of the material in the envelope falls back towards the binary system, we suggest 
th at the readjustm ent phase in the case studied here can be even more dramatic than in the case studied by suggested 
by Ilvanoval (|201lf ). In the case studied here there are two possible scenarios. In the first scenario the falling back 
material will be added to the remnant material and will be included in the readjustment phase. This material will 
be added to the outer layer of the left over envelope, which has the potential to expand and swallow the companion. 
In the second scenario, where there is enough angular momentum, a circumbinary disk is formed and the binary 
separation decreases as discussed above. Overall , our new finding of the existence of a large amount of fallback 
material suppo rts the conclusion of Ilvanoval (|2011l ) of a higher percentage of immediate post-CE mergers, although 
Ivanoval (|201lf ) mechanism works for giants with non-degenerate core. 



We emphasize the role played by the circumbinary disk. At early stages of the CE the spiraling in is relatively 
slow and the rate of energy release by the spiraling-in companion is low. At that stage the convection in the AGB 
envelope efficiently transfers the energy outwards. Most (but not all) of the extra energy is radiated away at this 
early phase. In the final phase of the CE, when the companion is deep inside, energy released at a very high rate, 
under the assumption used here. In that case energy is carried by the expelled envelope. This stops the spiraling- 
in process. Here we suggest that further spiraling-in can take place due to interaction with a circumbinary disk. 
The circumbinary disk is an efficient mechanism to absorb the binary orbital angular momentum. In addition, 
the circumbinary disk dissipates the binary orbital energy and radiates it by expelling mass and radiation from 
its surface. This phase ends either with a merger or with a very small orbital separation. A final small orbital 
separation will be interpreted as a small «ce ■ 

The value of ac e will be small even if our assumption of a rapid release of energy at the end o f the CE 
is not satisfied, but r a ther the release of ene rgy is a long process (e.g., ~ 1 yr, ISandquist et al.lll998t ~ 10 yr, 
De Marco et al.1 120031 ) . iDe Marco et al.l (|2011l ) suggest that if the spiral- in of the companion takes longer than the 
dynamical time, as in that case of long time energy release, the thermal energy of the AGB will be used to help 
unbind the envelope. Therefore, they conclude it is expected that the companion will spiral-in inward to the radius 
predicted by taking «ce = 1 
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4.2. The SN la scenario from an early merger 



Here we discuss some possible outcomes of WD- WD (double degenerate) me rger, where the two component 
are m ade of CO. The WD- WD merger is one of the scenarios that lead to S Ne la (Iben fc T u tukovlll984i IWebbink 



200 



3) assume the merger 



1984J). Most simulations and calculations of double degenerate merger (e.g., lYoon et al. 
to occur a long time after the CE phase, when the two WDs are already cold. In the scenario we propose here, the 
merger occurs within the final stages of the CE, while the core is still hot. Below we speculate that some of these 
will become SNe la after a very long time dictated by angular momentum loss. The speculation is based on studies 
of double degenerate mergers that we describe below. We then discuss their application to the scenario proposed 
here. 

In analyzing the outcome of the double degenerate merger we must take into account the following, (i) In 
the merger process the more massive WD remains at the center of the merger product, while the lighter WD is 
destructed and becomes the envelope, and possibly a thick disk around the center, (ii) The post-AGB core is hot as 
it had no time to cool, while the companion WD is colder as it had the time to cool since its formation. 

In general, the merger of two WDs can end in one of three basic products (e.g., Yoon fc Langer 20051) : SN la, 
core collapse to a neutron star (NS), or a more massive WD. In the first two cases, where the mass of the merger 
is above the critical limit, the most important factor in determining whether there will be an explosion as a SN la 
or not is the nature of the carbon ignition. In the third case, the mass of the merger is below the critical limit. We 
now discuss the three cases. 

(1) In the first case, the companion mass is larger than the post-AGB core mass, and therefore the companion 
will become the cold core of the merger (assuming it had enough time to cool) and the destructed hot post-AGB 
core will be its envelope. We note that before the merger but already in the CE phase, the center of mass will be 
closer to the companion, and the envelope will arrange itself as if the companion w ere the giant core. In this case 
carbon ignition is more likely to occur during the merger processes and be off-center ( Saio fc Nomotolll985 ). Most of 



the WD merger models indicate that a non- violent off-center ignition that occurs during the merger process, burns 
carbon and oxygen into oxygen, neon, and magnesium. As not enough energy remains to unb ind the WD when the 



Saio & Nomotol 


1985; 


Mochkovitch & Livio 


1990; 


a NS will be formed ( 


e.g.. Yoon & Langei 


2005). 



case 



Yoon fc Langerll2005l) . It is inter esting to mention th e possibility that in such a scenario 



planets can form around the NS from SN fallback material ([Hansen et alJl2009| ). 



However, in many cases no carbon ignition will occur during the merger process (|Yoon et al.l 120071 ). and if the 
mass reaches the critical mass a SN la might occur after the WD slows down. 

(2) In the second case, the companion mass is smaller than the post-AGB core mass, and therefore the center 
of the product will be hot. 



Yoon et al.l (|20071 ) raised the possibility that a hot remnant core is less likely to ignite carbon off-axis (in the 
envelope) during WD- WD merger. The reason is that a hot WD is larger, such that its potential well is shallower 
and the peak temperature of the accreted destructed- WD material is lower. Hence, in such a case supercritical-mass 
remnant is more likely to ignite carbon in the center at a later time, leading to a SN la. The merge r remnant becomes 
a rapidly rotating massive WD, that can collapse only after it loses sufficient angular momentum. lYoon et al.l (|2007[ ) 
parameterized the angular momentum loss, and found that explosion can occurs as late as ~ 10 6 yr after the merging 
occurs. 

(3) The Third possibility is that the merger is below the critical mass. In this case the merger will remain a 
WD that will continue to cool, and neither a SN la nor collapse are expected. 

However, there is also an exceptional case. Ivan Kerkwiik et al. I (|2010l) suggest that in cases where the two CO 
WDs have an approximately equal mass below the critical limit, the WD merger will be fully mixed and hottest in 
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its center. Initially, the merger rotates differentially, but in later times some of the mass is expelled, and a uniformly 
rotating core near the mass-shedding limit is obtained. The structure is of a core which contains ~ 80 per-cent of 
th e mass, surrounded by a s ub-Keplerian, very dense, partially degeneracy- pressure supported disk. In the model 
of van Kerkwijk et al. ( 201ol ) the disk accretes at a high rate that leads to compressional heating that is likely to 
cause ce ntral carbon ignition. Th is ignition occurs at densities for which pure detonations lead to events similar to 
SNe la (jvan Kerkwiik et al.l 120101 ) . If such a scenario occurs in our case, then the explosion will be surrounded by 
the expelled common envelope. 



4.3. Core-degenerate merger at the termination of the common envelope 



Our new suggestion is that there are many double degenerate merger events where one of the components, the 
post-AGB core, is hot. The timescale for angular momentum loss might be very long (> 10 6 yr after merger) even 
when the remnant rotation is modest: not fast enough to slow down rapidly, but fast enough to influence the critical 
mass. We hence suggest that some SNe la might originate from a double degenerate merger that occurs at the last 
phase of the CE or shortly afterwards during the planetary nebula phase, rather than a long time after. We term 
this scenario core-degenerate (CD), as one of the WD has just emerged from the core of the AGB star. 

The immediate implication of processes that reduce the orbital separation at the final stages of the CE (and 
later form an early merger) is that the effective value of «ce is considerably smaller than 1. The reduction in the 
orbital separation can lead to merger immediately after the CE phase, or even before the entire envelope has been 
lost. In a some cases the final orbital separation will be very small, such that the merger will occur after the CE 
phase, but while the core is still hot. We now examine this case. 



As mentioned in section 14.21 lYoon et al.l (|2007f ) showed that if the merger has a hot remnant core it is less 
likely to ignite carbon off-center, and therefore a massive CO WD is formed. The reason is as follows. If the post- 
AGB core is still hot when it mergers with a lighter WD companion, the core radius is still relatively large (e.g., 
Koester &: Schoenberner 1986t Bloecker 1995 ). When the post-AGB core radius is larger, then the gravitational 
potential well on the surface, where the companion mass is accumulated, is smaller. Therefore, when the companion 
WD merges with the post-AGB core, the peak temperature of the merger will be in an outer radius, and in a lower 
temperature. So the hotter the post-AGB core, the colder the merger peak temperature will be, for a given core 
mass. An increase by a factor of 1.2 in the post-AGB core radius results in a factor of ~ 1.2 decrease in the mergers 
peak te mperature. Such a decrease might be enough to prevent off-center carbon ignition and the collapse of the 
merger (jYoon et al.ll2007j ). Instead, the merger product will later produce a SN la. 



We consider the condition for a core-degenerate merger while the core is still ~ 1.2 larger than its final radius. 
As a post-p ost-AGB core c ools its radius decreases to an asymptotic value of a cold WD. More massive cores cool 
faster (e.g., Bloecker 19951) . For a SN la progenitor we take a post-AGB core of M core ~ 0.7-0.8 M Q with a WD 



companion mass o f M2 < M rorB ~ 0.6-0.7 M . Such cores shrink to a radius of 1.2 times larger than the final radius 



in T C ool ~ 10 b yr ( Bloecker 19951 ). If by that time the companion WD will merge with the post-AGB core then it 
will be less likely to ignite carbon off-center during the merging process. 

There are two possible mechanisms to further reduce the orbital separation in cases where the circumbinary disk 
did not reduce the separation down to merger. These are emission of gravitational waves and energy dissipation by 
tidal interaction. The timescale for the energy dissipation by tidal interaction is much larger, even longer than the 
Hubble time (e.g., IWillems et alJl2010l ). We therefore consider emission of gravitational waves to be the dominant 
mechanism. 



In an eccentric orbit WD- WD bin ary, the reduction of the s emi-major axis due to emission of gravitational waves 
is faster, as the eccentricity increases ([Peters fc Mathewslll9631 ). The semi-major axis reduction and the eccentricity 
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increase are ( Peters fc Mathewsl 1963 ; see also equations (2) and (3) in Ignatiev et al. 2001 ) 



64 G 3 M cme M 2 (M c , 



M, 



c 5 a 3(! _ e 2)7/2 



73 
24* 



37 
96* 



and 



where c is the speed of light. 



304 G 3 M core M 2 (M cmc + M 2 )e 



15 



c 5 a 4 (l-e 2 ) 5 / 2 



121 

304 f 



(27) 



(28) 



For the core-degenerate scenario to occur, the semi-major axis reduction due to emission of gravitational waves 
should result in a merger, within a period t 2 < t coo \. The value of t 2 , which satisfies t 2 » t\ (where t% is the time 
period of the circumbinary disk - binary interaction) is very sensitive to the value of e from which the GW emission 
process starts. This value of e is determined by the final eccentricity in the circumbinary disk - binary interaction 
process, which depends on t\. Namely t 2 strongly depends on t\. Our results are presented in Fig. [7] 



We repeated our calculation (with x = 2; equtionJM]) for different values of qd,o, and found that for qd$ > 0.1 
the two components will merge within a short time after the CE (Fig. [7]). As we take disk wind into account, we 
get that the process of emission of gravitational waves is required for merging to occur for 0.1 < qd.o < 0.12, but 
due to the interaction with the disk it will occur in a considerably shorter time than it would have taken without 
it. This shorter time can fulfil the requirement that when the core-degenerate scenario will occur, (namely that the 
companion WD will merge with the post-AGB core while it is still hot) and there will be a central ignition that 
results in a SN la. We also repeated the calculation for % = 1, namely with the disk wind having a zero terminal 
velocity. In this case we find that higher circumbinary disk mass (larger qd,o) is required for the companion WD to 
merge with the post-AGB core while the later still hot. 

This means that stars with massive envelopes (M env > 5 M Q ) that engulf a W D companion , might end the CE 
phase with a very close WD- WD binary. This is in agreement with the results of Livio fc Riessl ( 2003 ). who found 
that for a merger to occur the AGB must be massive. The WDs might merge before the hot WD (that was the 
post-AGB core) had time to cool much. We note that as tow is short, the probability of observing such systems 
is low. But, observationally, there are indeed post-CE systems that have periods much shorter than simulations 
predict (|De Marco et al]|2008l 2009). Our result that qce < 1 can account for this discrepancy. 



It is in place to mention that recent studies also reach the co nclusion that acE < 1- iDe Marco et al.l (|201lf ) 
perform simulations for many sets of binary parameters (see also iPassv et al.l 120111 ) . They find that there is a 



possib le negative correlation between the mass ratio of the two stars and the value of qce- As IDe Marco et al 



([20111 ) include the thermal energy (which is half of the binding energy according to the virial theorem), as we do 
here, they also find that the values of a qv. are below 1 (in case s where the thermal energy is not included, the values 
of acE they obtain are larger than 1). IZorotovic et al.l ([20101 ) also come into the c onclusion that acv, < 1. Their 
statistical analysis of observed post-CE binaries give «ce = 0.2-0.3, but contrary to lDe Marco et al.1 (|201ll ) they do 
not find any correlation between acE and the mass of the companion. 



5. SUMMARY 

We conducted an analytical study of the ejection of a common envelope (CE) under some simplifying assump- 
tions. Despite the many 3D numerical hydrodynamical simulations of CE evolution, such an approach has its own 
merit, as it points to the basic physics behind the finding that some of the envelope remains bound to the system. 

Most of the binary gravitational energy in a CE is released in a short time during the final stages of the CE 
phase. Here we took this time to be much shorter than the dynamical time in the outer regions of the envelope. 
Therefore, the energy release process resembles a blast wave propagating from the center outwards. We discuss two 
types of blast waves (section [3]) , which practically lead to the same result. In the first type the energy release is 
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instantaneous, i.e. an explosion, while in the second type there is a continuous energy deposition over a short time. 
We used a self similar solution to track the blast wave in each type as it propagates through the AGB envelope. 
We found that part of the ejected envelope stays bound to the binary system, i.e., it does not reach the escape 
velocity. This material is expected to fall back towards the center. It is very likely that due to angular momentum 
conservation and further interaction of the fallback gas with the binary system, a circumbinary disk will be formed. 
The spherically symmetric geometry of the self similar solution does not allow the introduction of angular momentum 
to the solution of the blast wave. 

As discussed in section 14. 1[ the binary semi-major axis decreases rapidly during the short circumbinary disk 

- binary interaction period. During that time the power of the system greatly exceeds the Eddington limit. Most 
of the energy goes to blow the wind, but some goes to radiation. It is therefore possible that such systems will be 
observed by the transient increase in their luminosity. 

We showed that interaction of the binary system with the circumbinary disk will further reduce the orbital 
separation. Consequently, in the alpha-prescription, we find that effectively «ce < 1. In many cases a merging will 
occur immediately after the dynamical phase of the CE. In many cases where the companion is a WD, it will form 
with the remnant post- AGB core of the AGB (or RGB) core a WD- WD system with very small separation. In cases 
where the mass of the circumbinary disk is large enough, we suggest that a core-degenerate (CD) is likely to occur 

- an early merger of the WD and the post- AGB core will merge while the post- AGB core is still hot. Such a merger 
is more likely to explode later as a SN la if the mass is above critical and the AGB's core is more massive than the 
WD companion (see section!?]). The super-critical mass WD is stabili zed by the rapid rotati on. Only after it slows 
down it will explode as a SN la. We speculate, based on the results of lYoon fc Langerl (|2005l ). that this time can be 
as long as > few x 10 9 yr, and that the core-degenerate scenario is another channel leading to a SN la. 



Our finding that effectively «ce < 1 (see also Ilvanovall201 II ) . can explain the recent findings of iDe Marco et al 



([201 II ). IDe Marco et al.l ([201 if ) find that the value of oce they deduce from observations is much smaller than 
what their numerical simulations of the CE phase give ( De Marco et al. 20081 2009, 2011). We strongly encourage 
numerical study of the CE phase to follow the gas and its angular momentum after ejection in order to find the 
amount of gas that falls back , and whether a circumbinary disk is formed. 
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Fig. 4. — The solution of an intense explosion blast wave, for an instantaneous release of energy inside the AGB envelope. We use AGB 
envelope mass M env = 4.2 Mq and a stellar radius of i?* = 310 Rq. The companion mass is M2 = 0.6 Mq. The physical quantities 
are plotted at the time the reaches the stellar surface (t en d = 9.47 days, upper panel), and at half of that time (t = i cn( j/2 = 4.73 days, 
upper panel). The plotted physical quantities are the density of the envelope p, the velocity v, the sound velocity v B , the escape velocity 
v CBC , the pressure p, the mass M inner to a radius r (including the core mass), and the total energy dE in a shell dr. In the upper panel 
we find that the intersection of the post shocked velocity v with the escape velocity v CBC occurs at a radius of ~ 150 R© (marked with a 
small arrow). The mass inner to this point does not reach the escape velocity and falls back to the binary system. We find the fall- back 
mass to be M lb ~ 0.22 Mq. 
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Fig. 5. — Physical properties as a function of the companion mass M2: the spiraling-in timescale at the final orbital separation 
(rj), the energy deposited in the envelope Eo (equation I13I I, the constant a (equation |12t . the fall-back mass Mr,, and the blast- wave 
propagation time to envelope expel time ratio t Bn d/Tf ■ 
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Fig. 6. — The solution of a blast wave with continuous energy release from t = to i = Tf ~ 0.95 days. A blast shell is formed, starting 
from the end of the void and ending at the Shockwave. For explanations to the legend see the caption in Fig. [4] 
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Fig. 7. — The variation of the semi-major axis a (in solar units) and eccentricity e as a function of time, in the core-degenerate 
scenario. We take \ = 1 m equation J25}, and an initial eccentricity eo = 0. The variation occurs in two phases. In the first phase the 
circumbinary disk - binary interaction reduces a and causes an increase in e. The final values of a and e in this phase are the initial 
values in the second phase - semi-major axis reduction and eccentricity increase due to emission of gravitational waves. This phase lasts 
a time interval T2, which strongly depends on the value of e at the beginning of this phase, and hence on t\. The first phase terminates 
when one of two things occur: the circumbinary disk is completely depleted, or when the binary system merge (occurs when e ~ 1). 
We find that the first phase is enough to create a merger for q^fi > 0.12. If the first phase has not resulted in a merger (because the 
circumbinary disk is not massive enough to result in a strong interaction with the binary), then the second phase of GW emission can 
reduce a to the extent of merging. We find that for 0.1 < qdfi < 0.12 the second phase will result in a merger within T2 < r coo i ~ 10 yr. 
According to the core-degenerate scenario, such an early merger that will explode as a SN la. For qd < 0.1 a late merger is formed, and 
such a merger is not expected to explode as SNe la. 




Fig. 8. — The same as Fig. but taking jf = 1 in equation {25}. Namely, with the disk wind having zero trminal velocity. In this 
case we find that higher circumbinary disk mass (larger q^ o) is required for the binary stars to undergo a fast merge. 



